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ABSTRACT 

We present a spectroscopic survey of the giant stellar stream found in the halo of 
the Andromeda galaxy. Taken with the DEIMOS multi-object spectrograph on the 
Keck2 telescope, these data display a narrow velocity dispersion of 11 ± 3 km s -1 , with 
a steady radial velocity gradient of 245 km s -1 over the 125 kpc radial extent of the 
stream studied so far. This implies that the Andromeda galaxy possesses a substantial 
dark matter halo. We fit the orbit of the stream in different galaxy potential models. 
In a simple model with a composite bulge, disk and halo, where the halo follows a 
"universal" profile that is compressed by the formation of the baryonic components, 
we find that the kinematics of the stream require a total mass inside 125 kpc of A/125 = 
7.5±i'j x 10 n M o , or M 125 > 5.4 x 10 11 M at the 99% confidence level. This is the 
first galaxy in which it has been possible to measure the halo mass distribution by such 
direct dynamical means over such a large distance range. The resulting orbit shows 
that if M32 or NGC 205 are connected with the stream, they must either trail or lag the 
densest region of the stream by more than 100 kpc. Furthermore, according to the best- 
fit orbit, the stream passes very close to M31, causing its demise as a coherent structure 
and producing a fan of stars that will pollute the inner halo, thereby confusing efforts 
to measure the properties of genuine halo populations. Our data show that several 
recently identified planetary nebulae, which have been proposed as evidence for the 
existence of a new companion of M31, are likely members of the Andromeda Stream. 



1 INTRODUCTION 



Stellar streams represent the visible remnants of the merging 
process by which the halos of galaxies are built up. By study- 
ing these streams we can attempt to unravel the formation 
of galactic halos, seeing when, how and how many small 
galaxies arrived and were incorporated into large galaxie s 
<Helmi fc White! Il999t IJohnston. Sackett fc Bullock! 1200 it) . 
Streams are also of great interest as probes of the large- 
scale mass distribution of the dark halos they reside i n 
Jjohnston et alJll99St llbata et al]l2001al; IZhao et alJll99Sft . 
This utility stems from the fact that streams from low-mass 
disrupting stellar systems trace the orbit of their progeni- 
tor, giving a means to constrain the tangential motion of the 
stars in the stream. The stars must move along the stream, 
and the magnitude of the tangential velocity must be such 
that, when the star is integrated along the orbit, it ends 
up with the same velocity as the stars that are currently 
downstream on the same orbit. 

Giant stella r streams may well be common structure s 
around galaxies JPohlen et alJl2003t iMalin fc Hadlevlll997ft . 

Indeed, the most conspicuous feature in the halo of the 
Milky Way is the giant rosette stream originating from 



the Sagittarius dwarf galaxy, which contains approximately 
half of the high latitude (|6| > 30 °) intermediate ag e stars 
at distances greater than 15 kpc lllbata et alJl2001al . liool : 
iMaiewski et alJl2003lh It would appear that the Milky Way 
has not incorporated into the Halo a more massive galaxy 
than the Sagittarius dwarf over the last ~ 7 Gyr. 

This naturally leads to the question of whether the 
MW has unusual feeding habits. To answer this, we 
have undertaken a large photometric study of M31, us- 
ing the wide field cameras at the INT and CFHT tele- 
scopes to resolve stars over the ent ir e disk and inner halo 
of that galaxy lllbata et alJ l2001bl: iFerguson et alJ 120021 ; 
iMcConnachie et alJl2003t) . This has given us an unprece- 
dented panoramic view of the large scale and small scale 
structure of a disk galaxy. The analysis of this huge dataset 
is still in progress, but it has already yielded some surpris- 
ing results regarding the incidence of substructure in the 
halo of M31. The most prominent of these substructures is a 
stream-like over-den sity of stars near the minor axis of M31 
lllbata et al.ll2001blh at first sight a facsimile of the Sagit- 
tarius Stream around the Milky Way. The red giant branch 
(RGB) stellar density in the halo increases on average by 
a factor of two in the on-stream regions and is statistically 
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Figure 1. The left hand panel shows the locations of the observed fields on the plane of the sky, in standard coordinates (£,»;). Asterisk 
symbols denote field centres of, from South to North, Fields 1-8 and 12-14 of our photometric survey of the Andromeda Stream with 
the CFHT12K camera; while the large rectangles display the size of the CFHT12K camera pointings. The Keck2 DEIMOS targets 
are shown as small d ots in the centre of Fields 1, 2, 6 and 8. Triangles represent the positions of the planetary nebulae identified by 
iMorrison et alj l2003t) . The inner ellipse demarcates the approximate limit of the visible disk of M31 at 2° (= 27kpc radius), and the 
outer ellipse shows a segment of a 50kpc radius ellipse flattened to c/a = 0.6, corresponding to the approximate limit of our INT 
survey. The right hand panel presents the radial velocity measurements (full circles), as a function of distance from M31 along the 
declination direction [= 780 tan(»;) kpc]. The positions of M32 and NGC 205 are also shown. The black line shows a straight-line fit 
[ v h{ii) = —4244.8 tan(?7) — 610.9kms~~ to the high negative velocity edge of this diagram. Stars that are more than 100 kms" 1 away 
from this line are marked with open circles. The two large squares show the chosen points for our "back-of-the-envelope" calculation, 
detailed in Section 3. The dashed, dot-dashed and dotte d lines show the orbital paths in simple Kepler, logarithmic and NFW potentials, 
respectively. The velocities of the IMorrison et alJ 12003ft PNe are represented again with triangles. 



significant at the 50-sigma level. Interestingly, this stream 
points toward the Andromeda satellites M32 and NGC205, 
and is aligned with the outer isophotes of NGC205, sug- 
gesting a relationship between the Andromeda Stream and 
these two dwarf galaxies. If this interpretation is correct, 
the stream has to be the result of previous interactions with 
M31, as there is otherwise not enough time to spatially sep- 
arate it from either of the two dwarf galaxies. 

The proximity of M31 provides us with an opportu- 
nity to undertake a spectroscopic survey of individual stars 
within the stream; indeed Andromeda offers the only ex- 
tragalactic giant galaxy in which such a study can be un- 
dertaken with current instrumentation, as in more distant 
systems such halo substructure would be smeared into very 
low surface brightness features. 

One useful property of the stream is that it is on a 
highly radial orbit. It passes very close to the centre of M31, 
where a comparison of the RGB tip of the stream with 
that of M31 itself shows them to lie at the same distance 
jMcConnachie et alj20"03T) ; whereas in the furthest field that 
it has been detected to date, the peak of its giant branch is 
0.27 mag fainter, indicating that it is 106 ± 20 kpc behind 



M31 jMcConnachie et aljEooll . This fortuitous alignment 
close to the line of sight, allows us to measure directly the 
potential gradient over > 100 kpc, and hence measure the 
halo mass. 

The layout of this paper is as follows: section 2 presents 
the spectroscopic survey of the Andromeda Stream, the re- 
sults of which are used in section 3 to constrain the mass 
of the dark matter halo of M31, and finally in section 4 
we discuss the significance of the results and present the 
conclusions of this study. Throughou t this work, we as- 
sume a distance of 780 kpc to M31 dStanek fc Garnavichl 
1998) and a systemic ra dial velocity of —300 kms 
ide Vaucouleurs et aljll99ll) . 



2 THE ANDROMEDA HALO SURVEY 

To understand the nature of the substructures detected in 
our panoramic halo surveys of M31, we undertook a follow- 
up programme to measure the kinematics of individual RGB 
stars with the DEIMOS spectrograph on the Keck2 10m 
telescope. On Sep 29-30 2002, we used the instrument in 



The course of the Andromeda stream 3 
line of sight depth (kpc) 




2 750 800 850 900 



f (degrees) Heliocentric distance (kpc) 

Figure 2. For convenience, the left-hand panel reproduces the chart shown previously in Figure 1. The right-hand panel shows a sideways 
view of the stream, drawn to the same scale as the left-hand panel, which displays the line of sight depth of the fields, together with the 
positions of M32 and NGC 205. The disk of M31 is highly inclined to our line of sight (12°5); the thick line is a schematic representation 
of a disk of radius 27 kpc inclined at 12°5. Evidently the stream orbits close to the plane of the Andromeda galaxy. 



the high-resolution configuration with the 12001/mm grat- 

o o 

ing (giving access to the spectral region 6400 A to 9000 A) 
to obtain high quality spectra of 768 stars in 9 fields in 
M31. This has given an order of magnitude improvement in 
statisti cs over previous radial v e locity surveys of the M31 
system llleitzel fc Guhathakurtal J2002J). The present contri- 
bution examines data from four fields along the giant stellar 
stream. 

Red giant branch stars were selected for observation 
from our CFHT12K camera survey, which covers the re- 
gion of the sky shown schematically on the left-hand panel 
of Figure 1. Stars were selected by choosing point-sources 
with f-band magnitudes between 20.5 < I < 22.0 and colours 
1.0 < V — I < 4.0. Both metal-poor and metal-rich popu- 
lations will be present in this selection region (see e.g. 
iMcConnachie et alJ 120031. The number of target stars per 
16.9' x 5' DEIMOS field is approximately 100. The slitlet 
widths were milled at 0"7 to match the median seeing, and 
the slitlet-lengths were always larger than 5", giving access 
to a good estimation of the local sky. 

The spectroscopic images were processed and combined 
using the pipeline software developed by the "DEEP" con- 
sortium. This software debiasses, performs a flat-field, ex- 
tracts, wavelength-calibrates and sky-subtracts the spectra. 
After extracting the spectra with a boxcar algorithm, the 
radial velocities of the stars were measured with respect to 
spectra of standard stars observed during the observing run. 
By fitting the peak of the cross-correlation function with a 
Gaussian, an estimate of the radial velocity accuracy was 
obtained for each radial velocity measurement. The accu- 
racy of these data are astonishing for such faint stars, with 
typical uncertainties of 5kms -1 to lOkms" 1 . Finally, M31- 
centred radial velocities are calculated by adding 300 km s -1 
to the Heliocentric values. 

In the 4 stream fields 184 stars were measured with 
radial velocity uncertainties less than 25kms _1 . The right- 



hand panel of Figure 1 shows the subset of 125 stars that 
have heliocentric velocity Vh < — lOOkrns -1 , as a function 
of projected distance dg(= 780 tan [77] kpc) along the decli- 
nation direction. For clarity, stars with Vh > — lOOkrns , 
which are primarily Galactic, have been omitted from this 
diagram. These four fields correspond to CFHT Fields 1, 
2, 4 and 8, and their locations on the sky are marked in 
the left-hand panel. The Heliocentric distances of the fields, 
as measured from the magnitude of the Tip of th e RGB of 
the Stream population iMcConnachie et alJliooll . are dis- 
played in the right-hand panel of Figure 2. The outermost 
field lies at (£ = — 2°0, r\ = — 4°1) at a Heliocentric distance 
of 886±20 kpc, that is, it extends out to 125T17 kpc from the 
centre of M31. This is the first time that it has been possi- 
ble to measure the kinematics of a stellar stream over such a 
huge range in galactocentric distance. An immediately strik- 
ing aspect of these data is the sharp edge to the distribution 
at negative velocities in Figure 1, with a velocity gradient 
that appears to be almost a straight line. The velocity dis- 
tribution near this edge is very narrow, as demonstrated by 
Figure 3, where we show the distribution of velocity offsets 
from the straight line in Figure 1. We compare the data 
between — 40kms _1 to 40kms _1 in Figure 3 to a Gaus- 
sian model using the maximum likelihood technique, and 
find that the dispersion is 11 ± 3kms _ . The distribution 
appears slightly skewed to positive velocities; defining skew- 

(see e.g..lPress et alJll992Tl. we 

find that the distribution of 82 stars in Figure 3 that have 
|w| < lOOkrns -1 have a skewness of 0.51. For comparison, 
the standard deviation of the skewness of samples of 82 stars 
drawn from a Gaussian distribution is 0.43, so the observed 
sample is not significantly skewed. Some of the stars in the 
skewed tail may be contaminants from the halo of M31, and 
others may be stream stars more distant along the line of 
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Figure 3. The distribution of velocities about the straight line 
displayed on the right-hand panel of Figure 1. The stream has 
a narrowly-peaked velocity distribution at the position of that 
straight line, with dispersion 11 ± 3kms _1 . 



sight. The nature of these objects may become clearer when 
N-body simulations are fit to the stream. 



3 CONSTRAINTS ON THE DARK MATTER 
HALO 

Before fitting an orbit model in a realistic galaxy poten- 
tial through the position and velocity data, it is worth in- 
vestigating what can be learned from a simple "back-of-the 
envelope" analytic calculation. We assume that the global 
potential at the distance of the stream is approximated 
by a spherical potential, and that the stream follows the 
orbit of the centre of mass so all stars have the same to- 
tal energy E = ^(r) ± \v(r) , where v is the three di- 
mensional velocity of a star at radial position r. If we as- 
sume further that the orbit is radial, which is a reason- 
able approximation given the distance information in Fig- 
ure 2, then v(r) = v(/3)/ cos(7), where v is the observed 
(one-dimensional) radial velocity, as a function of angu- 
lar distance j3 = cos -1 [cos(£ ) cos(?7)] = sin _1 [7 



L780kpc Sin (7)] 

along the stream, and 7 is the projection angle onto the 
line of sight. The M31-Sun-Field 1 triangle has an angle of 
cos _1 [cos(-2°0) cos(-4°l)], and lengths 780 kpc (M31-Sun), 
886±20kpc (Sun-Field 1), and 125 ±17 kpc (Field 1-M31). 
Thus the cosine of the projection angle onto the line of sight 
is 005(7) = 0.868 ± 0.040. We take 2 data points (£r = 
2°0,r ? i = -4°1) and (£ 2 = 0°7, r? 2 = -1°6), for which the 
de-projected Andromeda-centric distance and de-projected 
velocity are: (ri = 125 ± 17 kpc, «i = — 9±13kms~ 1 ), (r 2 = 
rising)/ sin(r?i) = 49±18kpc,u 2 = -222 ± 13kms _1 ), 
to represent the stream. The uncertainty on r\ is calcu- 
lated from the projection and distance uncertainties, while 
the uncertainties on Vi and « 2 are estimated from the 
11 ± 3 kms -1 dispersion combined with the projection un- 
certainty. These two points correspond to the extremities of 
Fields 1 and 6 (we do not choose Field 8, as this region is 
much deeper in the potential, where the disk contribution is 
significant). In calculating these de-projected velocities, we 
have assumed that the tangential velocity of M31 is zero. 



It has long been suspected that this tangential velocity is 
small given that there is no other large galaxy in the Lo- 
cal Group to p rovide significant torqu e to the Milky Way 
- M31 system jKahn fc Woltierlll959l) . iLvnden-Bell fc Linl 
dl977t) have also pointed out that this tangential motion 
must be small, as there is otherwise not enough time for 
M31 and the Milky Way to have almost completed an or- 
bit about each other in the age of the Universe, as is re- 
quired by the Local Group tim ing arguement. The result of 
lEinasto fc Lvnden-Belll (I1982T) is consistent with this pos- 
sibility; they find that the transverse motion of M31 is 
60 ± 30kms _1 , under the assumption that M31 and the 
Milky Way have equal and opposite angular momenta. To 
provide a crude assessment of the effect of the possible trans- 
verse motion of M31 on the parameters derived from our 
dataset, we also consider in the analysis below the conse- 
quence of a transverse motion of amplitude 300 km s -1 (i.e. 
equal to the radial velocity component) in the direction par- 
allel to the Stream. 

For a Kepler potential, the two chosen distance and 
velocity data points imply a central mass of 4.6 ±0.7 ±0.3 x 
10 11 Mq (the second uncertainty quantifies the effect of a 
300 kms -1 transverse velocity). A slightly more realistic case 
is a logarithmic potential $ = in 2 log(r 2 ± r 2 ), which has 
a circular velocity that asymptotes to v c at large radius r. 
Following iDehnen fc Binnevl 1^998) , we adopt a core radius 
of r c — 3 kpc for this model. In this case, 



(«2 



«i 2 )/ln( 



r 2 c +rj 
ri + r\ 



which gives a circular velocity of v c = 162 ± 15 ± 8 kms . 
We also investigated a simple model in which M31 is 
an NFW potential * = -GM S log (1 ± r/r s )/r, where 
M s is the mass within 5.3r s , and r 3 is the scale radius. 
Assuming the characteristic relation between halo mass 
and concentration in a z — LCDM cosm o logy c = 
15.0 - 3.3 logfM^nn /10 12 M^h' 1 ) iBullock et~ai] (l200ll) . we 
fmdM s = 4.4±1.2±0.3x 10 11 M , (M 200 = 8.0±2.1±0.5x 
lO 1:L M ,r 2 oo = 194±16±4kpc, v 2m = I36± 11 ± 3 kms" 1 , 
r s = 13.1 ± 1.4 ± 0.4 kpc). For comparison, the mass inside 
125 kpc is 7.6±1.2±0.3 x 10 11 M Q for the logarithmic model 
and 6.4 ± 1.3 ± 0.3 x 10 11 M Q for the NFW model, respec- 
tively. The projected distance-velocity behaviour of these 
three toy models is compared to the velocity measurements 
in Figure 1. It is no surprise that the Kepler potential over- 
predicts the radial velocity close to M31, but both the simple 
logarithmic and NFW halo models manage to approximate 
the stream velocity profile well. 

To constrain more realistic models of the poten- 
tial, we adopted a recen t comp osite galaxy model by 
iKlypin. Zhao fc Somervillel ([2002) , and calculate the po- 
tential due to a sum of disk, bulge and halo mass dis- 
tributions. The bulge in this model has a mass of 1.9 x 
1O 1O M , while the disk has a mass of 7 x 1O 1O M , and 
a s cale length of 5.7 kpc. Fixing the disk and bulge with 
the lKlypin. Zhao fc Somervillel (I2002T) parameters (with the 
further assumpti on that the disk scale height is 400 pc, 
as sug gested by iGouldl Il994h . we investigate the spher- 
ical halo models compatible w ith our data. Following 
IKlypin. Zhao fc Somerviilel {2002), we take into account the 
adiabatic contraction of the halo due to the settling of 
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Figure 4. A multidimensional view of the two best-fitting orbits in the more realistic galaxy model of lKlypin. Zhao fc Somervilld £2002). 
The kinematics in Figure 1, combined with the distance information of Figure 2, require that the stream orbit is currently moving in 
towards M31, gaining velocity from Field 1 to Field 6 as the stream approaches M31. The full line shows the orbit projection on the sky 
for the best-fitting orbit, in the best fitting potential, when only the data in Fields 1—8 are taken into account; in this case the total mass 
of M31 inside 125 kpc is M125 = 7.5^ 2, g X 10 11 Mq. The dashed line however, shows the best fit orbit in the best fitting potential when 
we also account for the distance data in Fields 12 and 13. Now the mass of the model is M125 = 1.5 ±0.1 x 10 12 Mq inside 125 kpc. In 
both cases the curves display the centre-of-mass orbits, integrated for ±0.75 Gyr from the position of the centre of Field 6; the arrowheads 
show the direction of motion. Panel (a) shows the location on the sky of the two orbits superimposed on the chart previously shown in 
Figure 1; panels (b) and (d) present projections of the distance-velocity relation (symbols as in Figure 1); and finally panel (c) shows a 
sideways view of the orbits (as in Figure 2). 



baryons into the disk and bulge, which alters the potential 
of the inner galaxy significantly. 

We consider a ha lo of a given mass, with the appro- 
priate concentrati on l|B ullock et alJ l2 00ll). and com press 
the halo using the iKlvpin. Zhao fc Somervilld <2002l) rela- 
tions. We then calculate the resulting potential by multipole 
expansion. Using an "AMOEBA" minimization algorithm 
JPress et al.lll992l) we launch orbits in the potential, itera- 
tively improving on guessed values for the initial position 
and velocity of a test particle on the orbit. A \ 2 goodness of 
fit statistic is minimized between the observed sky position, 
line of sight depth and radial velocity of the stars in our 
survey and the values predicted by the orbit model. This of 
course assumes that the stream follows the orbit of a test 
particle, an assumption that is reasonable if its progenitor 
was a relatively low-mass dwarf galaxy, which is consistent 
with the measured 11 ± 3kms _1 velocity dispersion in the 
stream. 

We assume that the orbit must pass through the field 
centres on the sky to within 0.15° (2 kpc); this corresponds 
to our estimate of the uncertainty in the measurement of 



the location of the central peak along the stream. To con- 
st rain the distance of the stre am, we take the measurements 
of M cConnachie et alJ (|2003) together with their 20 kpc dis- 
tance uncertainties (these uncertainties include an estimate 
of the systematic error in the distance measurement). The 
fit is also constrained with the measured radial velocities; 
we take all data-points within 100 kms -1 of the straight- 
line fit in Figure 1 (those that are shown with filled circles) , 
and adopt the fitted velocity dispersion of llkms" 1 as the 
expected uncertainty in the velocity fit. 

The relative likelihood of halo mass models can be ana- 
lyzed by comparing the likelihood of the orbit fit as a func- 
tion of halo mass. However, we find that the result depends 
strongly on whether the distance data of Fields 12 and 13 
are included in the the fit (the stream RGB population was 
not detected in Field 14). 

We first investigate the consequences of rejecting 
the distance data i n Fie lds 12 and 13. In those fields 
iMcConnachie et alJ l)2003|) may have detected the metal- 
rich thick disk, or a warp in the disk instead of the actual 
stream itself. Without confirmation from radial velocities, 
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Figure 5. The paths, in M31-centric coordinates, of the two 
best-fitting orbits previously shown in Figure 4. In this co- 
ordinate system, the Sun is located at (xM3i>yM3li z M3l) = 
(—762, 0, —169) kpc. The upper panel shows a face-on view of M31 
(the limits of a 27 kpc radius disk is indicted with a thick circle), 
while the lower panel is an edge-on view of the plane of M31. 
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Figure 6. The circular velocity curve of the two best-fitting po- 
tentials. The high-mass model (dashed line), based on the as- 
sumption that the stream continues to Fields 12 and 13, greatly 
over-predicts the observed rotation curve between 20-30 kpc. The 
M31 rotation curve data are reproduced here with sol id dots (from 
the compilation bv lKlvpin. Zhao fc Somervillell2002fl . 



circular velocity of the potential is shown with a dashed line 
in Figure 6. This orbit does not manage to fit the large line 
of sight distances in Fields 1 and 2 as well as the previous 
case (although the discrepancy is not statistically significant 
given the large error bars on these distance data) . The main 
shortcoming is that the circular velocity of this model at 
galactocentric distances below 30 kpc is at odds with the 
data displayed in Figure 6. The circular velocity of the halo 
alone (without the bulge and disk) is higher than the data 
allow. This model is therefore not acceptable. 



we cannot be certain that the stream continues into that 
area of M31. Fitting only the data in Fields 1 to 8, we find 
that the most likely model has a total mass inside 125 kpc of 
A/125 = 7.5 x 10 11 M Q ; the likelihood drops by a factor of e _1 
(i.e. la) for M 125 = 1.0 x 10 12 M Q . We also reject at the 99% 
confidence level a mass lower than A/125 = 5.4 x 10 11 Mq. 
The most likely orbit in this model galaxy potential is shown 
as a full line in the projections displayed in Figures 4 and 
5, and its circular velocity curve is displayed in the same 
way in Figure 6. This orbit manages to fit the distance in- 
formation in Fields 1-7 very well, but it is peculiar in being 
exceedingly radial (peri- to apocentre ratio of 70). Indeed, 
the orbit has a tangential velocity in Field 1, relative to 
M31, of 5kms~ , that is, the orbit has only a very small 
non-radial velocity component. This may be partly respon- 
sible for the good agreement of the present mass estimate 
with those derived above for the three toy models in which 
radial orbits were imposed. In the present analysis we have 
again assumed that the tangential motion of M31 is negligi- 
ble. However, the effect of a 300 km s^ 1 tangential velocity 
along the direction of the Stream alters the derived mass 
inside 125 kpc by only 0.5 x 10 11 M . 

If we decide to also fit the distance data in Fields 12 
and 13, the most likely model has a much higher total mass 
inside 125 kpc, with A/125 = 1.5±0.1±0.25 x 10 12 M Q (again, 
the the second uncertainty value reflects a tangential motion 
of M31 of 300 kms" 1 ). The most likely orbit in this potential 
model is shown with a dashed line in Figures 4 and 5, and the 



4 DISCUSSION AND CONCLUSIONS 

In this first kinematic study of the giant stellar stream in 
the Andromeda galaxy, we have been able to measure the ra- 
dial velocity gradient along the stream from the outermost 
field currently probed, 125 kpc distant from the centre of 
M31, down to an inner field ~ 20 kpc from the centre of that 
galaxy. Over this huge distance, the (projected) radial veloc- 
ity changes by 245 km s" 1 , implying a de-projected velocity 
difference of ~ 280 km s -1 . This velocity gradient is used to 
obtain a zeroth order analytic estimate of the mass of the 
halo, which for simple halo-only galaxy models such as a log- 
arithmic halo or an NFW halo, implies a mass inside 125 kpc 
of A/125 = 7.6±1.2xlO n M Q and A/125 = 6.4±T3xlO n M , 
respectively. In both cases, the uncertainty associated with 
a possible tangential velocity of M31 of 300 km s -1 , is ~ 5%. 

We also investigate more realistic solutions, allowing 
the stream to have a non-radial orbit, and taking a galaxy 
model that is the sum of a dis k, bulge and dark halo 
flKlvpfn. Zhao fc Somervillel l2002T l . The dark halo of this 
model is a perturbation on a spherical NFW model, to ac- 
count for the adiabatic contraction of the dark matter as 
the baryonic components form. If we disregard the distance 
data in Fields 12 and 13, since we cannot be certain that the 
stream is present in those regions, the most likely mass of 
this galaxy model is A/125 = 7.5^ 2 '3 x 10 11 Mq, with a lower 
limit of A/125 = 5.4 x 10 11 M Q (at 99% confidence). This 
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Figure 7. The spreading of orbits after a close encounter with 
M31. The thick-line orbit reproduces the best-fit orbit in the re- 
alistic galaxy potential previously shown in Figure 4, integrated 
from the centre of Field 6. The thin-lines show similar orbits, 
starting from the same spatial position, but with velocities per- 
turbed by a random offset drawn from a Gaussian distribution of 
dispersion llkms - 1 . This fanning-out of stars on nearby orbits 
will lead to the disappearance of the Andromeda Stream, and 
implies that it is a transient phenome non. The sma l l squa re on 
the minor axis shows the position of the lBrown et all J2003ll ACS 
field (not drawn to scale!). All other markings are as in Figure 1. 



result is fully consistent with the zeroth order analytic esti- 
mates discussed above, and suggests that the derived mass 
is not very sensitive to the adopted mass model. Further- 
more, there is a reasonable agreement of the resulting rota- 
tion curve with the kinematics of previously-observed disk 
tracers (see Figure 6), despite the fact that the model was 
only fit to the kinematics of the stream. This agreement is 
due in part to the fact that we have taken previously-fitted 
models for the disk and bulge, but the halo contribution to 
the total rotation curve dominates beyond ~ 20kpc, and 
it is at these large distances that we have fitted the model 
to the stream data. This confers further confidence on the 
derived mass model. The uncertainty on the mass estimate 
due to the possible tangential velocity of M31 is likely not 
very large, approximately 7% for a tangential velocity of 
300 km s" 1 . One of the main uncertainties in the present 
analysis is the flattening of the halo, which we have not ex- 
plored, as the current data set does not provide sufficient 
constraints. Future studies, fitting N-body simulations to a 
larger kinematic sample of stream stars, can be expected 
to improve the mass estimate and also constrain the halo 
shape. 

It is only recently that measurements of the mass of 



M3 1 beyond the edge of its gaseous d isk have been possi- 
ble. ICourteau fc van den Bergbl ill 9991) analysed the veloci- 
ties of 7 Andromeda satellites, finding a total mass of 13.3 ± 
1.8 x 10 11 Mq. In contrast, using a larger sample of 10 satel- 
lite galaxies, 17 globular clusters and 9 pl anetary nebulae 
as test particles, lEvans fc Wilkinson! (120001) found that the 
most likely total mass of M31 is 12.3lg 8 x 10 11 M , approxi- 
mately half of their Milky Way estimate of 19tfj x 10 11 M . 
With improved radial velocities of the M31 satellites they 
were later able to reduce their M31 mass uncertainties, 
findi ng a value of the total mass of ~ 7.0L°3 5 5 x 10 11 M 
lEvans et alj EoOO). which is fully consistent with our re- 
sult of M125 = 7.5± 11 x 10 11 M Q (defined within a radius 
of 125 kpc). However, a definitive statement of the relative 
masses of M31 and the Milky Way awaits an improved mea- 
surement for the Milky Way. 

The best fit orbit in the best fit potential is prograde 
and in the region where it is currently observed, it lies close 
to the plane of M31. Thus it appears that the orbit of the 
Andromeda stream is peculiar in being extremely radial, 
passing very close (within 2 kpc) of the centre of M31. This 
requires very special initial conditions. The stream stars, 
which are spatially narrowly confined in Fields 1 to 8, will 
diverge dramatically upon passing close to M31 to form a 
low-density fan-like structure, since orbits that deviate only 
slightly from the orbit displayed in Figure 4 on the plunging 
part of their course will take very different paths after being 
accelerated around the centre of M31, as demonstrated in 
Figure 7. The fanning-out of the stream is likely to confuse 
efforts to measure the metallicity an d age of the M31 ha lo; 
for instance, the recent discovery bv lBrown et alJ l2003t) of 
a young halo component in M31 from main-sequence fitting 
of an extremely deep ACS field may be due to stream con- 
tamination (see Figure 7). 

Thus the stream may be in the process of vanishing as 
a coherent structure, providing a supply of metal-rich stars 
into the halo. This also suggests that the stream was re- 
moved from its progenitor less than an orbital period ago 
(the pericentre to pericentre period of the continuous-line 
orbit in Figure 4 is 1.8 Gyr), as we would otherwise not ob- 
serve the structure as a stream. The ephemeral nature of the 
stream implies that the progenitor must have survived until 
~ 1.8 Gyr ago. As we discuss further below, the progenitor 
was probably of low mass, implying that the rate of decay 
of its orbit due to dynamical friction was slow, so it followed 
(or continues to follow) an orbit close to the current orbit of 
the stream. However, any dwarf galaxy on the derived orbit 
must have experienced extreme tides as it repeatedly passed 
close to the centre of M31. One option is that the progenitor 
was a very dense dwarf galaxy that was sufficiently robust to 
survive the huge tides. This brings to mind M32 as a can- 
didate, though detailed numerical modeling is required to 
examine this possibility. The alternative option is that the 
progenitor of the stream deflected off another halo object, 
sending it plu nging i nto th e current orbit, analogous to the 
suggestion by IZhaol ll998l) to explain the longevity of the 
Galactic satellite Sagittarius. 

However, the connection with M32 presents some diffi- 
culties. Although M32 appears to reside in the stream, its 
velocity is markedly different. For M32 to be associated with 
the stream would require it to be at a different phase in the 
orbit (either lagging or trailing). Furthermore, the low ve- 
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locity dispersion of the stream would appear to preclude 
M32 as its pr ogenitor (which has a v ~ SOkms" 1 outside 
of the nucleus. Ivan der Marel et, al.lll994h . though this can- 
not be confirmed without a detailed dynamical study. The 
case for association of NGC 205 with the stream also ap- 
pears weaker given these kinematic measurements, since the 
best fit orbit does not overlap with it in phase-space. Future 
studies may allow us to examine this issue in more detail 
by following the stream beyond the region currently probed 
with kinematics. The velocities of the stream stars presented 
here also shed light on the recent ide ntification of a possibl e 
new companion to M31, And VIII iMorrison et al.l [2003 ) . 
The radial velocities of those planetary nebulae, which are 
located close to M32 (see Figure 1), appear to have radial 
velocities consistent with an extrapolation of the Stream, as 
they lie close to the straight line in the right-hand panel of 
Figure 1. This would suggest that And VIII is most likely 
part of the Andromeda Stream, situated in the re gion of 
highest over-density reported bv llbata et al ] fcOOlbh . How- 
eve r, it is interesting to note that the radial velocities of 
the iMorrison et all i2003t) planetary nebula sample tend to 
increase towards the North (showing an apparent positive 
gradient in the right-hand panel of Figure 1), whereas the 
Stream stars have a negative velocity gradient. The connec- 
tion between the two structures therefore merits to be exam- 
ine d more careful l y. Du ring the refereeing process, a study 
by iMerrett et alJ i2003l) was presented which also investi- 
gates the planetary nebulae around M31. However, though 
their survey has also detected PNe in a region at the base 
of the stream (near our Field 8), their interpretation is in- 
consistent with the kinematics of the stars reported here. 
The direction of motion of the orbit they derive is oppo- 
site to ours, and the path of their orbit, which intercepts 
many M31-disk PNe is substantially different to the orbit 
that we have fitted. It is possible that their finding reveals 
the presence of another kinematic structure in Andromeda. 

The majority of the stream stars that were surveyed 
have a narrow velocity dispersion of 11 ± 3kms _1 , though 
slightly skewed to positive velocities. The fact that our line of 
sight looks down the stream (see Figure 2), so that we prob- 
ably see stars over a range of distance along the line of sight, 
and hence at different phases in the orbit, will tend to ren- 
der the observed velocity dispersion higher than the intrinsic 
velocity dispersion. This indicates that the progenitor was 
most likely a low mass dwarf galaxy. The Milky Way satel- 
lite Sagittarius, wh ich has a velocity dispersion of 11 kms -1 
fabata et al.lll997l) . also has a gigantic stellar stream, but 
with a larger velocity dispersion of 20kms _1 iYannv et alJ 
2003). However, it is unclear at present whether the lower 
velocity dispersion measured in the Andromeda stream com- 
pared to the Sagittarius stream implies that its progenitor 
was of lower mass than Sagittarius, or not. 
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